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Abstract. Since the Hipparcos mission and recent large scale surveys in the optical and the near-infrared, new 
constraints have been obtained on the structure and evolution history of the Milky Way. The population synthesis 
approach is a useful tool to interpret such data sets and to test scenarios of evolution of the Galaxy. We present here 
new constraints on evolution parameters obtained from the Besancon model of population synthesis and analysis 
of optical and near-infrared star counts. The Galactic potential is computed self-consistently, in agreement with 
Hipparcos results and the observed rotation curve. Constraints are posed on the outer bulge structure, the warped 
and flared disc, the thick disc and the spheroid populations. The model is tuned to produce reliable predictions in 
the visible and the near-infrared in wide photometric bands from U to K. Finally, we describe applications such as 
photometric and astrometric simulations and a new classification tool based on a Bayesian probability estimator, 
which could be used in the framework of Virtual Observatories. As examples, samples of simulated star counts at 
different wavelengths and directions are also given. 



Key words. Galaxy : stellar content 
- Galaxy : structure 



Galaxy : general - Galaxy : evolution - Galaxy : kinematics and dynamics 



1. Introduction 

The population synthesis approach aims at assembling 
current scenarios of galaxy formation and evolution, the- 
ories of stellar formation and evolution, models of stellar 
, atmospheres and dynamical constraints, in order to make 
^ ' a consistent picture explaining currently available obser- 
^ - vations of different types (photometry, astrometry, spec- 
troscopy) at different wavelengths. 

The validity of any Galactic model is always ques- 
tionable, as it describes a smooth Galaxy, while inhomo- 
geneities exist, either in the disc or the halo. The issue is 
not to make a perfect model that reproduces the known 
Galaxy at any scale. Rather one aims to produce a useful 
tool to compute the probable stellar content of large data 
sets and therefore to test the usefulness of such data to 
answer a given question in relation to Galactic structure 
and evolution. Modeling is also an effective way to test 
alternative scenarios of galaxy formation and evolution. 

The originality of the Besangon model, as compared to 
a few other population synthesis models presently avail- 
able for the Galaxy, is the dynamical self-consistency. The 
Boltzmann equation allows the scale height of an isother- 
mal and relaxed population to be constrained by its veloc- 
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ity dispersion and the Galactic potential. The use of this 
dynamical constraint avoids a set of free parameters quite 
difficult to determine: the scale height of the thin disc at 
different ages. It gives the model an improved physical 
credibility. 

In section 2, we show how the model describes the 
Galactic populations, as they are accounted for in the 
model and how constraints on galaxy evolution are in- 
corporated. In section 3 we describe current and future 
applications of such a model and give a sample of com- 
parisons between model predictions and existing data. In 
section 4 we discuss ongoing plans and future improve- 
ments. 



2. Galactic populations and Galaxy evolution 

The scenario of evolution of the Galaxy slightly evolves 
as new data become available (new wavelengths, deeper 
samples, more accurate observations, ...), and as physical 
processes are better understood. Model predictions have 
to be checked and eventually model parameters have to 
be refined. One should always take care that the model 
stays self-consistent and that the overall agreement is con- 
served. In this section we describe the overall inputs in 
an up-to-date model of the Galaxy, how accurately the 
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physical parameters are constrained, and which data are 
used to constrain them. This results in a description of the 
Galactic populations consistent with the current scenario 
of Galaxy evolution. 

The main scheme of the model is to reproduce the 
stellar content of the Galaxy, using some physical assump- 
tions and a scenario of formation and evolution. We essen- 
tially assume that stars belong to four main populations 
: the thin disc (sect. 2.1), the thick disc (sec. 2. 2), the stel- 
lar halo (or spheroid, sect. 2.3), and the outer bulge (sect. 
2.4). White dwarfs are taken into account separately but 
self-consistently (sect. 2.5). A model of distribution of the 
extinction is also assumed (sect. 2.6). 

The modeling of each population is based on a set of 
evolutionary tracks, assumptions on density distributions, 
constrained either by dynamical considerations or by em- 
pirical data, and guided by a scenario of formation and 
evolution. 

In iBienavme et al.l l|l987albf l we have shown for the 
first time an evolutionary model of the Galaxy where the 
dynamical self-consistency was taken into account to con- 
strain tlie_disc_scale height and the local dark matter. 
Then iHavwood et alJ l)l997f) used improved evolutionary 
tracks and remote star counts to constrain the initial mass 
function (IMF) and star formation rate (SFR) of the disc 
population. The thick disc formation scenario has been 
studied using photometric and astrometric star counts in 
many directions, which also provided its velocity ellipsoid, 
local density, scale height, and mean mctallicity. These 
physical constraints led to a demonstration that the prob- 
able scenario of formatio n was by an accre tion event at 
early ages of the Galaxy l)Robin et all Il996|) . 

More recently, the Hipparcos mission and large scale 
surveys in the optical and the near-infrared have led to 
new physical constraints improving our knowledge of the 
overall structure and evolution of the Galaxy. These new 
constraints, now included in the present version of the 
model, are presented below. 



2.1. Thin disc 

In the scheme of computing the population synthesis 
model, the stellar content at each epoch has been com- 
puted from standard parameters such as the initial mass 
function (IMF), a star fo rmation rate (SFR). a nd a set 
of evolutionary tracks (see iHavwood et al.L Il997l and ref- 
erences therein). The disc population has been assumed 
to evolve during 10 Gyr, age derived from the white 
dwarf lumino sity function with an accuracy of about 15% 
ijWood k Oswalt, .1998.) . Sets of IMF slopes and SFRs 
have been tentatively assumed and tested against star 
counts. The preliminary tuning of the disc evolution pa- 
rameters against relevant observa t ional data was exten- 
sively described in IHavwood et alJ l)l997|) . We present be- 
low further improvements in the disc description from new 
observational constraints. 



2.1.1. Hipparcos inputs and the dynamical 
self-consistency 

Thanks to the Hipparcos mission, new and more accu- 
rate measurement s have been mad e of the local den- 
sity and potential llCreze et al.[ ll99?^. luminosity function 
jJahreiss & WiclcnI. 11997^ and kinematics ( Gomez ct al] . 
fl997h . These new constraints have been included in the 
modeling as follows. 

First, the new Hipparcos luminosity function is used 
to constrain the IMF at low masses in the solar neigh- 
bourhood. The IMF is modeled by a two-slope power law 
(table 1). At high masses, the slope a = 3 is slightly 
higher than the Salp eter assumption, as determined by 
IHavwood et al ] |ll997l) . At the low mass end the slope has 
been r evised to a=1.6 (instead of 1.7 in Havwood et al.j 
to give a better fit to the Jahreifinew luminos- 
ity function from the revised Catalogue of Nearby Stars 
(2001, private communication). The new IMF slope gives 
a better fit to the local luminosity function in the range 
6 < My < 13 as shown in figure 1. At the lower lumi- 
nosity end, the local sample is still too poor to ascertain 
the IMF slope. This slope is in good agreement wit h inde- 
pendent and detailed analysis from lKroupal l)2000|) which 
is in favor of a slope a = 1.3 ± 0.5 for the mass range 
0.08 to 0.5 Mq. Kroupa adopted four different slopes in 
four mass ranges. Contrary to us, in the mass range 0.5 to 
1.0 Mq, Kroupa adopts a = 2.3 ± 0.3, that is a Salpeter 
IMF. We estimate that presently available field data do 
not give enough constraints to adjust an IMF with four 
different slopes depending on the mass range. We have 
rather chosen to rely on a simpler modeling (two slopes) 
of the IMF for the thin disc population, until complemen- 
tary constraints are available. 

Second, as Hipparcos has provided new constraints on 
the local kinematics, the age-velocity dispersion relation 
from Gomez et al. (1992) is introduced (cf sect. 2.1.4). 

Third, the local density of interstellar matter has been 
revised. From iDamel ( 1993fl the local surface density of 
HI-I-H2 is about 6 Mq pc~^ with a considerable uncer- 
tainty of about a factor of 2. With a scale height estimated 
to 140 pc, the local mass density is 0.021 M0 pc~'^. This 
is about half the local density used in the previous version 
of our model. We notice however that the uncertainty is 
still very large. 

Adding up the mass density from stars, interstellar 
matter and local contribution of the dark halo (see ta- 
ble 2), we end up with a total local mass density of 0.0759 
Mq pc~^, in agreement wit h the dynam i cal va lue 0.076 
±O.O15M0pc-3 obtained by ICreze et all l)l998|) from a 
detailed analysis of the force on a sample of A and F 
stars observed by Hipparcos. 

The evolution model gives the present distribution of 
stars in a column of unit volume centered on the sun as a 
function of intrinsic parameters (age, mass, effective tem- 
perature, gravity, metallicity). Since the evolution model 
does not yet account for orbit evolution, we redistribute 
the stars in the reference volume over the z axis according 
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Fig. 1. Luminosity function in the V band in the solar neighbourhood. Crosses with Poisson errors: Jahreiss determi- 
nation from the revised Catalogue of Nearby Stars; Lines: model LF assuming an IMF slope of a=1.7 (solid line) and 
1.6 (dashed line) at m<lM0. 



Table 1. Age, metallicity ([-^l) (mean and dispersion about the mean), radial metallicity gradient (dex/kpc), initial 
mass function (IMF), and star formation rate (SFR) of the stellar components. 





Age (Gyr) 


[^] (dex) 


dR 


IMF 


SFR 


Disc 


0- 0.15 
0.15-1 

1- 2 


0.01 ± 0.12 
0.03 ± 0.12 
0.03 ± 0.10 




dn/dm oc m~" 






2-3 


0.01 ± 0.11 


-0.07 


a = 1.6 for m < IM© 


constant 




3-5 


-0.07 ± 0.18 




a = 3.0 for m > IM© 






5-7 


-0.14 ± 0.17 










7-10 


-0.37 ± 0.20 








Thick disc 


11 


-0.78 ± 0.30 


0.00 


dn/dm oc m"" '' 


one burst 


Stellar halo 


14 


-1.78 ± 0.50 


0.00 


dn/dm oc m""'^ 


one burst 


Bulge 


10 


0.00 ± 0.40 


0.00 


dn/dm oc m"^ -^" for m > O.7M0 


one burst 



to the age/thickness relation deduced from the Boltzmann rotation curve constr ains the shape and loca l density of 
equation. the dark matter halo ijBienavme et al.l Il987aj) . 



With these new inputs, the dynamical self-consistency 
process is applied as first described in iBienavme et alJ 
assuming that the Galactic potential is given by 
the sum of the stellar populations (those in the present 
model), and adding up the interstellar matter and the 
dark matter halo, we compute the potential through the 
Poisson equation. The disc is not mo deled wi t h a d ouble 
exponential but with Einasto laws ijEinastd. Il979(l . en- 
suring continuity and derivability in the Galactic plane. 
The formulae are given in table 3. The thin disc popula- 
tion is sliced into seven isothermal populations of different 
ages, from to 10 Gyr. Each sub-population (except the 
youngest one, which cannot be considered relaxed) has 
its velocity dispersion imposed by the age/velocity dis- 
persion relation. Then using equation ^ we compute its 
scale height (or axis ratio in case of ellipsoids), while the 



p{0) 



exp 



'w 



(1) 



A new total mass density is then computed. The pro- 
cess is iterated until the values of the potential and the 
disc thicknesses change by less than 1%. The resulting disc 
axis ratios e are given in table 2. The rotation curve of the 
best fit model is shown in figure 2. 

The local density of the dark matter halo amounts to 
9.9 X 10"^ MfT^pc"'^, about 20 % higher than the value 
determined by IBienavme et al while the mass 

density of our disc is lower. The difference comes from the 
fact that they included an "unseen mass disc" of local den- 
sity 0.01 MfT^ pc"'^ of which exi stence has been definitely 
ruled out bv lCreze et al.l l)l998f) . The dark halo serves to 
compensate for the lack of force coming from the disc, in 
order to adjust the rotation curve. This is the major mas- 
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Table 3. Density laws and associated parameter of the stellar components, — R"^ + ^ where R is the galactocentric 
distance, z is the height above the Galactic plane, and e is the axis ratio. Values of e are given in table 2. The disc 
density law is given here without the warp and flare. The corrections for these structures are given in the text (sect. 
2.1.3). do are normalization constants. For the bulge, x,y and z are in the bulge reference frame and values of N, 
xo,yo, zo and Rc are given in the text (sect. 2.4). 



density law 



Disc 



Spheroid 
Bulge 

ISM 

Dark halo 



po/do X {exp{-{a/hR^) ) - exp(-(a//iH_ ) )} 



with = 5000 pc, hfi_ 
po/do X {exp{-{0.5^ +a^/hj 



= 3000 pc 

J^/2)-exp(-(0.52+aV/i^ 



0/2 



)} 



with ftfl , — 2530 pc, hR_ — 1320 pc 



Thick disc po / do x exp ( 



po X exp (- 
with hii 
po/do X if-) 



R-R, 



hf 



)X(1 



cxp(3:i/h^) / \z\ N 

2500 pc, hz = 800 pc 



N X exp(-0.5 X rf) 

„2 



N X exp(— 0.5 X r^) x exp(— 0.5( 

with =yp^)TTiiWT{^^ 

po X exp{~^-^) X ea;p(-j^) 

with hB = 4500 pc, = 140 pc 
(i.+jl"/ Ra)'^) with -Re = 2697 pc and pc 



0.1079 



if age < 0.15 Gyr 
if age > 0.15 Gyr 

if \z\ < xi, a;i=400 pc 
if \z\ > xi 



if a < flc 
if a > Be 



flc = 500 pc 



y/x'^ + y2 < 
y/x^ + y^> Rc 



Table 2. Local mass density po of the stellar compo- 
nents, the dark matter halo and the interstellar medium 
(ISM). IF- velocity dispersion aw used for the dynamical 
self-consistency, and disc axis ratios e resulting from this 
process. The white dwarf mass density is computed as- 
suming a white dwarf mass of O.GMq. 





Age 
(Gyr) 


(M0PC-3) 


aw 

(km s~^) 


e 


Disc 


0-0.15 


4.0x10"-' 


6 


0.0140 




0.15-1 


7.9x10"^ 


8 


0.0268 




1-2 


6.2x10"^ 


10 


0.0375 




2-3 


4.0x10"^ 


13.2 


0.0551 




3-5 


5.8x10"^ 


15.8 


0.0696 




5-7 


4.9x10"^ 


17.4 


0.0785 




7-10 


6.6x10"^ 


17.5 


0.0791 




WD 


3.96x10"^ 






Thick disc 


11 


1.34x10"^ 








WD 


3.04x10"* 






Stellar halo 


14 


9.32x10"*^ 




0.76 


Dark matter halo 




9.9x10"^ 




1. 


ISM 




2.1x10"^ 







M0, excluding the central black hole and clusters. The 
mass of the interstellar matter amounts to 4.95 x 10^ 
MqIu our modeling. The stellar thin disc reaches a mass 
of 2.15X IQi" Mq, the thick disc 3.91 x 10^ Mgand the 
stellar halo 2.64x 10* M© (excluding eventual ancient halo 
white dwarfs which belong to the "dark halo"). The dark 
halo has a mass of 4.53 x 10^^ Moinside the sphere of 
50 kpc radius. The overall Galaxy mass to R<50 kpc 
amounts to 5.04x 10" Mq, and 9.97x 10" Mgwithin 
100 kpc. These values are in good agre ement with the 
detail ed analysis of the Galactic mass bv lSakamoto et all 
(2003), from the latest kinematics of Galactic satellites, 
halo globular clusters and horizontal branch stars, which 
estimates the Galactic mass to 5.5j]Q'4 x 10^^ Mq within 
50 kpc and about 10 x 10" Mq within 100 kpc. 

The density laws resulting from self-consistent adjust- 
ments of the Galactic potential are finally used to compute 
the stellar content all over the Galaxy, assuming that the 
disc IMF and SFR at the solar neighbourhood are suitable 
for the entire disc population. The inferred dark matter 
local density is also used to compute the contribution of 
potential ancient white dwarfs in star counts (see sect. 
2.5.3.) 



sive component at galactocentric distances larger than 5 
kpc. The question of the composition of this dark compo- 
nent is still open, probably mainly non baryonic, but not 
exclusively. We shall consider in sect. 2.5 the possibility of 
this component to be made of ancient white dwarfs. 

The overall model mass density ends up with the fol- 
lowing values: the outer bulge mass, coming from old pop- 
ulations in the central part of the Galaxy, is 2.03 x 10^" 



2.1.2. Atmosphere models and colour estimates 

The evolutionary model computes the distribution of stars 
in the space of intrinsic parameters: effective temperature, 
gravity, absolute magnitude, mass and age. This distri- 
bution transcribes the theoretical constraints coming from 
stellar and galactic evolution (mainly the IMF, the SFR, 
the evolutionary tracks, and the metallicity distributions) 
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At present, calibrations are available for the Johnson- 
Cousins UBVRIJHKL. The GAIA G magnitude can also 
be simulated. The Sloan SDSS and Stromgren photomet- 
ric systems are planned. 



10 20 30 

Galactocentric distance (kpc) 



Fig. 2. Circular velocity of the self-consistent 
model. Crosses ar e obse rvational constraints from 
ICaldwell fc Ostriked ^8^. Contributions from main 
Galactic components are also given. 



into observable statistics, with the use of suitable atmo- 
sphere models. 

The metallicity is estimated from the age of each com- 
ponent, assuming an age/metallicity relation and allow- 
ing for a Gaussi an dispe r sion about the mean. For the 
thin disc, we use iTwarod l)l98Cl() age/metallicity distribu- 
tion (mean and dispersion about the mean, see tabic 1). 
An alt ernative new metallicity distribution from jjaywoocl 
more centered on the solar metallicity, can also 
be used. However important to understand the scheme 
of the Galaxy chemical evolution, the differences between 
Twarog and Haywood distributions have only marginal 
impact on star count predictions as far as broad photo- 
metric bands are used. Conversely, it is not possible to 
settle this issue with ordinary photometric star counts. 

The theoretical parameters are converted into 
colours in various photometric systems through stel- 
lar atmosphere models co rrected to fit empirical data 
ijLeieune et all ll997L IT99I (the BaSeL 2.2 data base). 
While some uncertainties still remain in the resulting 
colours for extreme spectral types, the BaSeL data sets 
allow an efficient and consistent conversion of theoreti- 
cal parameters (effective temperature, gravity and metal- 
licity) to colors in various photome tric systems. For low 
mass stars, synthetic colours from iBaraffe et al.l ((l998l) 
have been used in the near-infrared in place of Lejeune 
calibrations. Empirical densities, infrared luminosities and 
colours of giants of spectral types M7 and later are intro- 
duced following the observations of .Guaii elmo (1993) . The 
mid-infrared modeling of these late M giants and AGB 
stars may be inaccurate. It will be revised in the near 
future. 



2.1.3. External disc 

Several studies have shown that the edge of the 
disc is detect ed at a g alactocentric distance o f 
about 14 kpc l|R,obin et all Il992l: iR^iphv et all Il99fit) . 
Howev er when analy si ng D ENIS near-infrared survey 
data l|EDchtein et al 

1 l^997^ in the Galactic plane, 
we have shown that the warp and flare of the 
external disc appear c l early in near-infrared star 
counts l|Derriere fc R.obi For a long time the 

warp and flares have b een clearly seen in HI data 
llHeiiderson et J]. Il982t iBurton and_1^e Linkel Hekkertl 
Il986t iBurt on. 1988; Diplas and Savage, 199l\j 
lar clouds (>Grabelsky et al 



40 



in molecu- 



1987;.Wouterlout et al 

May et all ll997l) . from OB stars llMivamoto et af 



199C; 



1988; 



Reed, 1996) a nd flnally in COBE dat a (Porcel et all 11997 
iFreiidenreick 119981: IPrimmel fc Soergel .200il) . The exis- 



tence of such a warp in old stellar populations puts con- 
straints on the time scale of this structure. However there 
is evidence that the stars do not follow exa.ctly the same 
warp as the gas I Diorgovski and SosinLll989tlPorcel et all 
I1997I) . Possible origins of the warps are still under discus- 
sion but may be resolved by suitable N-body simulations 
l|Binnevl [T992'). Warps may have originated from interac- 
tions between the disc and (i) the dark halo (if angular 
momenta are not aligned), (ii) nearby satellite galaxies, 
such as the Sagittarius dwarf or the Magellanic clouds, (iii) 
infalling interga lactic g as. From the analysis of their an- 
gular momenta, lBailinl\2003.) argues that the Milky Way 
warp and th e Sagit tarius Dwarf Galaxy may be coupled. 
iGarcia-Ruid l)2002|) have studied by N-body simulations 
the possibility for the warp to be caused by the tides of 
the Magellanic Clouds. They find that neither orientation 
nor amplitud e of the warp can be rep r oduce d in this way. 
Alternatively iLopez-Corredoira et all l|2002() propose that 
the warp caused by intergalactic accretion flows onto the 
Milky Way disc. 

Despite the uncertainty on the warp's origin, in the 
modeling process we have started with warp and flare pa- 
rameter values close to the ones observed in the gas and 
molecular clouds. Then, simulations have been compared 
with near-infrared data from the DENIS survey in a set 
of relevant directions, and warp and flare parameters have 
been adjusted. The parameters are describe d below : 

Fo r the warp, a tilted ring model such as IPorcel et alJ 
1I1997I) is considered. When computing density laws for 
a warped disc, the galactocentric coordinates {R,9,z) at 
galactocentric radii larger that i?warp, are shifted perpen- 
dicular to the plane by a value 

^warp — COs(^ ^max) 

with 0„iax and Zc being respectively the direction in which 
the warp is maximum, and the maximum elevation of the 
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warp at a given radius. We adopt a linear increase of the 
warp amplitude Zc, such as 

— Twaip(-^ -Rwarp) 



Following iBurtonI l|l988f) . we assume that the Sun lies 
approximately on the line of nodes of the warp 
90°). We adopt an amplitude similar to lCvuk et alJ fit 
(7warp = 0.18). Preliminary results obtained from the 
analysis of a part of the DENIS survey ijPerriere fc Robinl 
I2OOI: Derriera.2001) indicate a best value for the starting 
galactocentric radius of the warp to be i?warp — 8.4 kpc, 
close to Rq. 



As in iGvuk et al.l l)l999(l . we model the flaring of the 
disc by increasing the scale height by a factor /cflaro, beyond 
a galactocentric radius i?flarc, with 

fcflaro = 1 + 7fiarc(^ - -Rflaic) 

The surface density in the flaring disc is the same as in 
the normal one, leading to correction factors in the density 
laws. 

The amplitude is taken from iGvuk et alJ lll999t) 

(iflare = 5-4 X 10~^ kpc""'^). The radius at which the disc 
starts flaring lies beyond the solar circle. However lDerrier j 
found that the minimum radius of the flare might 
depend on the longitude considered. He determined a pre- 
liminary average value of i?flarc = 9.5 kpc which may be 
reconsidered when near-infrared surveys such as DENIS 
and 2MASS will be completely analysed. 

An example of the warp and flare effects on star counts 
in the near infrared is shown in figure 3. Squares are 
DENIS K star counts on a strip crossing the plane at 
1=266.1° and dashed lines are model predictions in same 
directions. The peak of the model predictions without 
warp and flare (top panel) is displaced from the data peak. 
The middle panel shows predictions including the warp 
but no flare. The Galactic warped plane is well placed 
but the peak is too sharp. In the bottom panel the model 
including both features (warp and flare) correctly repro- 
duces the star counts. 



2.1.4. Disc kinematics 



The d isc velocity ellipsoid is taken from iGomez et al.l 
l|l997l) . deduced from Hipparcos data analysis. Velocity 
dispersions as a function of age and galactocentric gradi- 
ents are given in Table 4. The asymmetric drift is com- 
puted following the relation : 

v;d = ^{^^+2x^^+(i-§)+(i-^)} 



2VlSR ^ P „ „ 

where ^ is the density gradient. The Sun velocities 
are Uq = 10.3 km s-\ Vq = 6.3 km s-\ Wq = 5.9 km 
s""'^, and Vlsr = 226 km s^^. The later value is obtained 
after fitting the mass model to the observed rotation curve 
(see sect. 2.1.1). Asymmetric drift values computed at the 
solar position are given in table 4. 



2.2. Thick disc 

The thick disc formation scenario that drives the model 
characteristics is a formation by a merger event (or sev- 
eral merger events) at the beginning of the life of the thin 
disc. This scenario well explains the observed properties 
of this population : The accretion body heats the stel- 
lar population previously formed in the disc making a 
thicker population with larger velocity dispersions. Due 
to the epoch of this event, the thick disc abundances are 
intermediate between the stellar halo and the present thin 
disc. Chemical abundance measurements imply that the 
duration of star formatio n in the t hick disc cannot be 
larger than 1 Gyr ( Mashonkina fc Gehren . 2001) . The es- 
timated epoch of thic k disc formation in the liter ature 
ranges between 8 Gyr l*Ibukivama fc Arimotoi l2002|) and 
13 Gyr jPetting er _e_t_alj , 20illj f rom c hemical consider- 
ations, while lOiiillen fc^ Garne^tl ll2001^ find 9 Gyr from 
kinematics. We have finally adopted an age of 11 Gyr for 
this population. This is just slightly older than the thin 
disc, and younger than the stellar halo. The choice of the 
slightly different age (by 1 Gyr) would not change signifi- 
cantly wide band photometric star counts. Other observa- 
tional constraints come from kinematics : the velocity el- 
lipsoid and asymmetric drift have been measured by radial 
veloc ity or proper motion distributions. Oiha et al. (1996, 
Il999l) determined the velocity ellipsoid of the thick disc 
from an analysis of a photo-astrometric survey in three di- 
rections at intermediate latitudes (center, anticentre, and 
antirotation). They found the velocity dispersions to be 
(cr^/, (Ty , (Tvf) = (67, 51, 42) in km with an asymmet- 
ric drift of 53 km s~^. These v alues are in close agreement 
with a recent determination bv lSoubiran et sl\ lj2003t) who 
find a rotational lag of 51 km s~^, and a velocity ellipsoid 
(63±6, 39±4, 39±4) km s^^. The asymmetric drift is not 
found to have any vertical gradient inside this population. 
The eventual radial gradient is still undetermined and is 
assumed to be zero in modeling. The absence of chemi- 
cal and kinematical gradients and the large gap in prop- 
erties between the thick disc and the stellar halo (scale 
height, rotational velocity) rule out the hypothesis that 
the th ick disc has been formed during a dissipational col- 
lapse ( Maiewskil 1 19931) and rather favour a formation by 
a strong heating of a preliminary thin disc induced by a 
merging of a satell ite onto the Galaxy l(Quinn et alJ . ll993l: 
iRobin et allll996(l . 

In order to model this population, a single epoch of 
star formation (that is a period of formation negligible 
compared to the age of the Galaxy) is assumed. Oxygen 
e nhanced evol utionary tracks have been taken from 
Ber gbush fc VandenBere. (1992X The horizont al branch 
(HB) has been added following models from iDormanI 
()1992|) . Stars of effective temperature between 7410 K and 
6450 K are assumed to be RRLyrae. There is no blue hori- 
zontal branch, as seen on metal rich globular clusters. The 
position of the HB is assumed to be 3.54 visual magnitudes 
brighter than the main sequence turnoff. 
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Fig. 3. Star counts in K along a DENIS strip in counts per 0.8° bin in declination (this strip crosses the plane at 
1=266.1°). Data are in solid lines, models in dashed lines. Top panel: flat disc. Middle panel: warped disc. Bottom 
panel: warped and flaring disc. 

Table 4. Velocity dispersions {aij,av,crw)i asymmetric drift V"ad at the solar position (see section 2.1.4) and velocity 
dispersion gradient —^w^j where ajj is expressed in km s^^ and R in kpc. 





Age 




av 


aw 




dR 




(Gyr) 


(km s-i) 


(km s-i) 


(km s-i) 


(km s-i) 


Disc 


0-0.15 


16.7 


10.8 


6 


3.5 
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8 


3.1 
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17.6 


10 


5.8 
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7.3 
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43.1 
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17.4 


14.8 
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43.1 
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85 
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The adopted thick disc metallicity is —0.78 dex, in 
agreement with in situ spectroscopic determination ft'om 
Icilmor c et al. (1995) a.nd photometric star coun t determi- 
nations ijRobin et all Il996t iBuser et all Il999l) . The low 
metallicity tail of the thick disc is not taken into account, 
but an internal metallicity dispersion of 0.3 dex is in- 
cluded. 

The thick disc density law has often been assumed in 
the literature to be exponentially decreasing perpendicu- 
lar to the Galactic plane. Tests of various laws have been 
performed by comparing data with model predictions in 
various directions. At present time, it does not seem pos- 
sible to obtain detailed constraints on the true shape, ei- 



ther exponential or scch^, because the distances are not 
known with sufflcient accuracy. Moreover the contribution 
of short distance stars to the star counts is often negligi- 
ble due to the small volume where they are counted. We 
have assumed that the thick disc shape is a truncated ex- 
ponential along the z axis: at large distances the law is 
exponential, at short distances it is a parabola (see ta- 
ble 3). This option ensures the continuity and derivability 
of the density law in the plane, to ease the computation 
of the thick disc contribution to the potential. 

The density law parameters and IMF adopted are well 
constrained by the analysis of deep wide field and near- 
infrared star counts at high and intermediate Galactic lat- 
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itudes l|Revle fc RobinL 120014 - The best fit model has a 
scale length hn = 2500 ± 500 pc, a scale height = 
800 ± 50 pc with xi=AOO pc, and a local density po ~ 
1.03x10-3 stars pc'^ or 7.6x10"'' Mq pc'^ for My < 
8, corresponding to 6.8% of the thin disc local density. 
The IMF is dN/dm oc m~° -^, the slope be i ng also con- 
strained by large data sets l|Revle fc Robinl.l2001a(l . The 
total density (excluding white dwarfs) is then 2.83x10"'^ 
stars pc~3 or 1.34 10"^ M© pc"^. This value does not 
account for binarity. Thus the true IMF should slightly 
steepen. The binarity correction should play a role only in 
very deep star counts. Such star counts are few and they 
remain uncertain due to strong galaxy contamination. 



2.3. Stellar halo 

We here distinguish what we call the stellar halo, or 
spheroid, from the dark matter halo. The spheroid is dis- 
tinct from the bulge in the model construction. However 
they may be related in their formation scenario. The 
spheroid, or stellar halo, is essentially old and metal poor. 
The bulge is also quite old but its abundance distribu- 
tion is sensitively different, with a metallicity distribution 
closer to the one of the thick disc population. On the other 
hand, the stellar halo should be considered distinct from 
the dark matter halo, which is probably mainly non bary- 
onic, although a small part of it can be formed of stellar 
remnants. We shall consider this point in sect. 2.5.3. 

Mainly two possible scenarios of formation of the 
spheroid are competing. The first one is a dissipational col- 
lapse of pre-galactic gas, the second one being based on ac- 
cretion of smaller pieces. Streams have already been found, 
related to known accretion events (Sagittarius galaxy, 
...) but the proportion of spheroid stars which can orig- 
inate from such streams is still not known. Therefore 
our standard model assumes a homogeneous population 
of spheroid stars with a short period of star formation, 
which should be a good representation on large scales. 
Comparisons with true d ata sho uld th en, by contrast, help 
finding existing streams. iBerebush fc VandenBerd l)l992|) 
oxygen enhanced models are used assuming a Gaussian 
metallicity distribution of mean -1.78 and dispersion 0.5 
dex. The age of the spheroid is assumed to be 14 Gyr, a 
value which may be slightly high, as constrained by th e 
analysis of the WMAP experiment ( Benne tt" et al.L l2003() . 
However, a change by 0.5 Gyr would have no effects on 
star count predictions. 

The density law and IMF of the spheroid were con- 
strained using dee p star counts at h igh and medium 
Galactic latitudes ("Robi n et all kOOd) . This global ad- 
justment of the spheroid shape has led to the following 
values: the density law (see table 3) is characterized by a 
power law index n = 2.44, a flattening e = 0.76 and a local 
density po = 2.185x10"^ pc^^ or 0.932 10"^ M© pc'^ in- 
cluding all red dwarfs down to the hydrogen burning limit, 
but not the white dwarfs. This local density has been com- 
puted with an IMF dN/dm oc m~°-^, similar to the globu- 



lar clusters and to the thick disc population (sec discussion 
below). The total local density depends a lot on the as- 
sumed IMF at low masses. With an IMF of dN/dm oc m~^, 
the total local density would be' 8.19 x 10~^ stars pc~'^. 
Model predictions are affected by the IMF slope value only 
at magnitudes fainter than V= 22 at the Galactic pole (see 
figure 4 and sect. 4). 

The relative density of the spheroid to the disc is about 
0.06% when one selects stars at My < 8^^. Yet this value 
ignores possible old white dwarfs. The ancient white dwarf 
contribution will be considered separately below. 

2.4. Outer bulge and disc hole 

Bulge formation in the Galaxy is still an open ques- 
tion. Several scenarios compete to explain the observed 
struct ures. The proposed scenarios are : primordial col- 
lapse (jfi ggen et al l Il962j) . hierarchical galaxy formation 
ijKauffmann et all Il994l) fearlv merging), infall of satel- 
lite galaxies (late merging), and secular ev olution of discs 
(starbursts induced by bar instabilities). iBouwens et alJ 
l|1999() argue that all these mechanisms could have been 
at wor k over the history of the Universe in forming 
bulges. In akasato fc Nomotiol 1 20031) have performed three- 
dimensional hydrodynamical N-body simulations of the 
formation of the Galactic bulge, and they find that two 
groups of stars are present in the bulge, 60% of the stars 
originate in sub-galactic mergers, and 40% come from the 
inner region of the disc. The two groups are chemically 
distinct. Their model is able to reproduce the metallicity 
distribution and the velocity dispersion-metallicity rela- 
tion. 

The observational signatures of the bulge population 
are the following : the bulge is certainly triaxial. The 
asymmetry is seen from COBE/DIRBE near-infrared data 
: the bulge is brighter a t posit ive lor igitudes than a t ne. 
ative ones. iDwek et"all l|l995l) and iFreudenreictf (1998 
have fitted parametric models to the deredc lened DIRBE 
data, but ex cluding low latitude region s. iBinnev et alJ 
l)l997t) . then i Bissantz fc GerhardI (|2002!) have obtained 
non-parametric models from the same data. The ob- 
served metallicity distribution of K giants is character- 
ized by a wide distribution, [Fe/H] r anging f r om -1 .2 to 
-^0.9 dex llMcWiniam fc RichL Il994^ . iMinnitil ((l996l) and 



iTiede and Terndrunri|1999l) have shown a clear relation 
between the velocity dispersion and the metallicity but 
this result may be due to a mixture of populations (halo, 
bulge and maybe thick disc) in their fields at latitudes be- 
tween -6 and -8. The bulge is rotationally supported, the 
projected mean rotation o f the bulge being of the orcl er of 
53.3 ± 0.3 km s"' kpc"' l)Tiede and Terndrunl [l999|) . 



' We apol o gize for a mistake in the published paper 
iRobin et al.L l^porf) where the given local densities of the 
spheroid were wrong by a factor of 2 

^ the ratio computed with low mass stars included would be 
completely dominated by the IMF slope used, parameter which 
is not strongly constrained 
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We here concentrate on the outer bulge, excluding the 
very central part of the Galaxy (inner bulge at longitude 
and latitude smaller than 1 deg) not yet well observation- 
ally constrained. 

As a first approxi mation, we have st arted with the 
Dwek G2 density law ( Dwek et all Il99fj) and have fit- 
ted the parameters on near-infrared star counts from the 
DENIS survey. A detailed description of the process will 
be given elsewhere (Picaud et al., in preparation). Here 
we give the main results from this work. 

The disc scale length is adjusted at the same time as 
the bulge parameters. The data towards the inner Galaxy 
clearly show a coexistence of both populations. However, 
since the Einasto formula for the disc accounts for a possi- 
ble hole in the middle of the disc, we take the opportunity 
to constrain this parameter and check how the bulge and 
disc populations coexist or are replaced by each other as 
a function of Galactic radius. 

The density laws are given in table 3. The fitted pa- 
rameters for the disc are the scale length hji^ and the hole 
scale length hfi_ . For the bulge there are 8 fitted param- 
eters (following Dwek notation): 

— orientation angles: a (angle between the bulge major 
axis and the line perpendicular to the Sun - Galactic 
Center line ), /3 (tilt angle between the bulge plane and 
the Galactic plane) and 7 (roll angle around the bulge 
major axis) 

— scale-lengths xq on the major axis, yo and zq on minor 
axes 

— normalization (star density at the center) TV and cut- 
off radius Re- 
in order to perform the Monte-Carlo si mulations, the 

bulge luminosity function is taken from iBruzual et alJ 
^ll997^ at solar metallicity and age of 10 Gyr. Alternatively 
Padova isochrones may be used in the future. A Salpeter 
IMF has been assumed (see tabic 1). This IMF slope seems 
to be well constrained for high mass and intermediate 
mass stars. This is however not the case at low masses, for 
which very deep observations inside the bulge are needed 
but difficult to get because of the crowding. From deep 
photometry obtained with the Hubble Space Telescope 
iHoltzman et al.l l)l998j) found a power-law mass function 
with a slope -2.2 for masses higher than 0.7, though close 
to the Salpeter IMF. At lower masses they found indica- 
tions of a break in the mass function, as seen in other 
populations, occuring at about 0.5 to 0.7 Mq. However, 
due to incompleteness and uncertain binary corrections, 
they do not strongly constrain the value of the IMF slope 
at low masses. If we would follow the general scheme of 
other populations, we would have taken either the low 
mass slope of the disc (0.6) or that of the thick disc and 
stellar halo (-0.5). The constraints being still poor we have 
kept the one-slope Salpe ter IMF for the t ime b eing. 

The DENIS survey ijEpchtein et all Il997|) has pro- 
duced specific observations in the bulge well suited for this 
analysis. We have used catalogues in J (1.25 fim) and Ks 



(2.15 fim) bands from a set of "batches", specific obser- 
vations made in the directions of the bulge (Simon et al, 
in preparation). Contrary to the strips done in the main 
survey, the batches are shorter in declination and easier to 
handle. About 100 windows of low extinction of 15'xl5' be- 
tween longitudes [-8°, -|-12°] and latitude s [-4°; -|-4°] were 
chosen using the map of extinction from ISchultheis et alJ 
|l999). The completeness limits were carefully established 
in each band. Star counts in K and J-K were simulated 
and compared with the data between magnitude K 8.5 
and 11.5 after extinction was fitted independently in each 
window. 

The model goodness-of-fit is estimated through a max- 
imum likelihood test. In order to explore the 10 dimen- 
sional space of parameters, sets of parameters were chosen 
by Monte Carlo selections, first from a uniform distribu- 
tion, then from a Gaussian distribution around the max- 
ima of likelihood. Fifty independent fits were made to test 
the robustness of the results. 

The best fit parameters (mean of the 50 trials and 
dispersions about the mean) are given in table 5. 

Thus, we obtain a very oblate outer bulge (sometimes 
called a bar in the literature), almost in the Galactic 
Plane with an angle from the Sun - Galactic Center line of 
about 10 degrees, well in agre ement with the result from 
iLopez-Corredoira et alJ l|2000|) who found an angle of 12 
degrees from independent star c ount data. The axi s ratios 
are in good agreement with the lDwek et al.l ||199,'t() deter- 
mination but they found a slightly higher angle, although 
less constrained, between the major axis and the line of 
sight to the Galactic center (20° ± 10°) compared to ours 
(11.1° ±0.7°). The good accuracy we achieved comes from 
the use of star counts, more sensitive to lower luminosity 
stars covering a larger depth inside the bulge and with a 
higher density, rather than integrated luminosities, domi- 
nated by the brightest stars, which may be biased by the 
deprojection. 

The overall mass in this outer bulge is 2.03 ± 
0.26 X 10^° M0. It includes all populations, even white 
dwarfs which contribute to about 26%. It does not in- 
clude the inner bulge or the central black hol e. This value 
is slightly higher than the one determined bv lDwek et alJ 
(1995) from COBE data 1.3 x IO^^Mq, but the latter 
is photometrically determined and does not account for 
white dwarfs. 

The overall fit of the central region leads to include a 
noticeable hole in the thin disc, which is very common in 
barred galaxies. The resulting disc scale length is hji^ — 
2.53±0.11 kpc and the hole scale length = 1.32±0.14 
kpc. 

The bar is considered as a separate population from 
the bulge, extending fur t her a,w ay (up to 1=27° as seen by 
Looez-Corredoira et al. Noo")). It is being studied and is 
discussed elsewhere (,Picaud et al., .2003.) . 
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Table 5. Bulge fitted parameters. Values from lDwek et al.l l)l995l) are mean va lues obtained from the 2.2 fj,m DIRBE 
data for their G2 density law. a, /3 and 7 are orientation angles following the iDwek et al. notation. xo,yo, zq 

are the scale lengths in the three axes. Rr is the cutoff ra dius and N the central number density. The values of (3 and 
Rc are not adjusted in the case of the lDwek et al fitting process. 



Parameter 


Dwek et al. value 


Error 


Our fitting value 


Error 


Unit 


a 


68.5 


4.8 


78.9° 


0.7 





13 







3.5° 


0.2 





7 


90.7 


1.5 


91.3° 


1.1 





Xo 


2.03 


0.23 


1.59 


0.12 


kpc 


yo 


0.470 


0.010 


0.424 


0.015 


kpc 


Zo 


0.580 


0.010 


0.424 


0.015 


kpc 


Rc 


2.4 




2.54 


0.13 


kpc 


N 






13.70 


0.7 


star.pc 



2.5. White dwarfs 

White dwarfs (WD) are the latest stage of stellar evolu- 
tion. Their density is a function of the past star forma- 
tion history and the initial mass function. Their luminos- 
ity function depends on the assumed cooling tracks and 
physical properties. In the modeling we assume all white 
dwarfs to be DA a nd use evol utionary tracks and atmo- 
sp here mode l s fromlBereeron et al. (J995|j), complemented 
by IChabrieJ l|l999^ for the very cool end, applicable to 
the halo. The former have been extensively compared with 
observations and have been shown to be reliable. The lat- 
ter models still have to be checked, as very few ancient 
white dwarfs are known yet. However, they include the 
best physical constraints and are supposed to be good 
starting points for simulating ancient white dwarfs. 

2.5.1. Thin disc 

The luminosity function (LF) of disc white dwarfs is 
known from systematic searches in the solar neigh- 
bou rhood. A g o od p reliminary LF was obtained 
by iLiebert et alJ l)l988|) . More recent determinations 
give larger local d ensities due t o new detections of 
cool white dwarfs: iHolberg et all ll2002t) fo und 5.5 ± 
0.8 X 10-3 star pc"'^ while iLeggett et al.l l|l998l) and 
lR,uiz fc Berge ron^ JioOl') foimd a lower limit of 5.6 x lO"'^ 
star pc-3, and Kn ox et aL (1999) 4.16 x 10^^ star pc^^. 
It should be noted that these values depend a lot on the 
number of faintest white dwarfs detected at My > 16, 
which a re only 3, and the way the luminosity function is 
binned. IWood fc Oswaltl l| 19981) estimated that the accu- 
racy of the white dwarf density is about 50% on a sample 
of 50 objects (the number of white dwarfs known in the 
local sphere of 13 pc is 53). 

To avoid Poisson noise from observational da ta, we 
have a dopted a luminosity function based on the IWoodI 
l)l992t) model for a disc of 10 Gyr, with a constant star for- 
mation rate, coherent with our evolutionary model. This 
is well in agreement wi t h the bright DA luminosity func- 
tion from lLiebert et ahl l)l988|) and is slightly higher in the 
coolest part of the observed white dwarf luminosity func- 
tion, allowing for incompleteness in the local sample. The 



total local density is then 6.6 xlO ^ star pc ^ and the 
mass density 3.96 xlO'^ M0 pc^^. 

2.5.2. Thick disc 

The thick disc white dwarf population has been com- 
puted by Garcia-Berro (private communication and 
iGarcia-Berro et al.l l|l999l )) assuming a Salpeter IMF. 
Very few thick disc white dwarfs have yet been identi- 
fied in the field. Hence, we have tentatively normalized 
the theoretical LF by fit ting model predictions on the 
lOppe nhcimcr et alJ ( 20011) sample of high velocity white 
dwarfs. Rcvlc et al.l 1 2001b) have shown that most of the 
Oppenhcimcr WD belong to the thick disc populations, as 
suggested by their kinematical properties. It gives a den- 
sity of white dwarfs relative to the red dwarfs of 20% that 
is 5x10^'* stars pc^'^. 



2.5.3. Halo white dwarfs 

Due to the microlensing results towards the Magellanic 
Clouds, the hypothesis that a significant number of an- 
cient white dwarfs are present in the halo is seriously con- 
sidered. The local contribution of this faint population in 
the solar neighbourhood is detectable and several studies 
have been conducted in order to find them. Thus, this hy- 
pothetical population has been tentatively included in the 
model. 

The current estimation of machos density from 
rnicrolensing events towards the Mag ellanic Clouds 
( Lasserre et al.[ 120001 lAlcock et all I2OO0I) has an upper 
limit of about 2 0% of dark matter halo made of machos. 
iFreese fc GraffI (|l997l) have shown the improbability of 
the machos to be red or brown dwarfs, leaving the pos- 
sibility for ancient white dwarfs to be the lenses. These 
results have stimulated the search for the local counter- 
part of the halo baryonic dark matter in the form of 
white dwarfs. However, the results give a quite smaller 
local density : iGoldman et alJ l)2002l) have shown that the 
local density of h alo white dwarfs cannot e xceed 5% of 
the dark halo and lOppenheimer" et alJ 1I2OOII) give a lower 
limit of 3%, although we do not think, from their kine- 
matics, that all Oppenheimer objects are really halo white 
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dwarfs l)Revle et al.U2001b|) . iMended l)2002|) has recently 
found a population of old white dwarfs in a proper mo- 
tion survey near the globular cluster NGC6397. Among 
his 7 peculiar white dwarfs he claims to have a halo white 
dwarf and 6 thick disc white dwarfs. If he is right, then 
the number of thick disc white dwarfs is huge compared 
with the number generally assumed and they have (for 
4 among 6) a surprising rotational velocity higher than 
the LSR, while one expects a significantly smaller rota- 
tional velocity for the thick disc population. The identifi- 
cation of these stars as belonging to the thick disc is then 
problematic. One should consider rather the possibility 
of them being disc white dwarfs which have un dergone a 
kick o ff in a mass-transfer binary, as proposed bv lBereeronl 
1)2003(1 . Considering the Mendez's halo white dwarf candi- 
date, the halo white dwarf mass density inferred by this 
object would be about half the density of the dark halo 
(4.25 ± 2.8 M0pc~^) if this star is a halo object. Its proper 
motion is still compatible with thick disc kinematics. 

Another recent search for high proper inotion ob- 
jects in the VIRMOS survey l|Le Fevre et all Il998|) on 
0.94 square degree adding to a 0.16 square degree in the 
SA57 field obtained also at the CFHT produces no detec- 
tion of high proper motion white dwarf in selection limits, 
where a 100% white dwarf halo aged 14 Gyr would place 
between 12.6 and 20.4 detectable objects (depending on 
the adopted halo IMF). This sets at the 95% probability 
level the upper limit of the halo fraction in the form of 
halo white dwarfs between 24% and 15%(Creze et al, in 
preparation). The thick disc white dwarf population seen 
by Mendez has not been seen in this Virmos field. 

IChabrier et al.1 l)l996|) have reconsidered the hydrogen 
white dwarf cooling sequence and recomputed their model 
atmosphere and color-magnitude diagrams. In order to 
build their luminosity function and det ection probabil i- 
ties, they considered two different IMFs dChabrierlFrnflfll) . 
They are mainly constrained by the nucleosynthesis which 
places limits on the amount of high mass elements rejected 
in the interstellar medium at early ages and by the fact 
that none were yet observed in white dwarf systematic 
searches. Of course there are not yet any direct obser- 
vational constraints on the true IMF of these first stars. 
Hence we have alternatively used both IMF to compute 
dedicated simulations. The local normalization is assumed 
to be 2% as a low conservative value compatible with an- 
cient white dwarf searches. The mean white dwarf mass is 
about 0.8 Mofor Chabrier's IMFl, and 0.7 Mefor IMF2. 
The local number density is deduced from these numbers 
and the local density of the dark matter halo as deter- 
mined by the dynamical constraints (sect. 2.1.1). We end 
up with a local density of 2.5 x 10~* star pc^'^ for IMFl 
and 2.8 x 10"'' star pc^^ for IMF2. 

2.6. Interstellar extinction 

The interstellar matter is mostly concentrated in the 
Galactic plane. Hence the intervening material can be 



modeled, at medium and high latitudes, as a double ex- 
ponential or an Einasto law, similar to the thin disc of 
youngest stars. However, when comparing model simula- 
tions with data in the Galactic plane, this smooth model 
is representative of the mean extinction only, and does 
not account for small scale variations of the extinction. 
At very low latitudes, it is better to apply, on a case by 
case basis, alternative distributions determined from ob- 
served colour-magnitude diagram s or from independent 
studies. iDrimmel fc Spergell l|200lll have built a 3D extinc- 
tion model, in which the distribution of extinction clouds 
in a thin disc and spiral arms is forced with the far-infrared 
emission due to the dust, as obser ved by COBE, or as de- 
termined bv lSchlegel et al.1 l(l998l) . An overall analysis of 
the near infrared star counts from recent all-sky surveys 
should provide in the near future a more accurate extinc- 
tion distribution due to dust in the Galactic disc. 

The reddening in each photornetric b a nd is obtained 
from the extinction law of lMathi'3 l|l99(|l . ICardelli et alJ 
il989) have shown that their extinction law, very closed 
to the latter, is well reliable over the wavelength range 
0.125 /J,m < A < 3.5/iTO and is applicable to both diffuse 
and dense regions of the interstellar medium, except for 
very dense clouds. We adopt a value of the commonly used 
parameter Ry = A(V)/E(B-V) = 3.1. This value may be 
inappropriate in some particular lines of sight in highly 
extincted clouds, which we do not intend to model in de- 
tail. Moreover the effect of the choice of Ry is negligible 
at A > 0.7/zm. 



3. Model simulations and applications 

3.1. Simulations 

Using the inputs described above, we compute the num- 
ber of stars of a given age, type, effective temperature, 
and absolute magnitude, at any place in the Galaxy. 
The assumed sun position is at a galactocentric ra- 
dius of 8.5 kpc (the lAU va l ue) and at 15 pc above 
the plane llFreudenreichL 19981: iHammerslev et all Il995t 
IPrimmel fc SnergelL I2OOI ). From the expected density of 
a given type of stars, a series of Monte-Carlo drawings 
are performed to make a sample of stars, the number of 
which is compatible with the expectation and includes the 
Poisson noise. Thus, when producing a series of catalogues 
in the same theoretical conditions one obtains variations 
in the number of stars compatible with the Poisson statis- 
tics. 

Apparent magnitudes and colours are then computed 
star by star for a given set of filters with the above men- 
tioned model atmospheres, reddened as a function of in- 
tervening extinction on the line of sight. The photometric 
bands currently available are the Johnson-Cousins system 
(UBVRIJHKL) and the GAIA G magnitude. Predictions 
at 7 and 15 microns for ISO bands, 12 and 25 microns 
IRAS bands are being tested, as well as the Sloan pho- 
tometric system. Photometric errors are included on each 
photometric band independently, the error being a func- 
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tion of the magnitude (as usually observed). This function 
can be a third degree polynomial, or an exponential as in 
equation 2. 



A + exp{C X m — B) 



(2) 



where m is the magnitude in the band considered, and 
A,B,C are fitted to the set of data with which the model 
predictions are going to be compared. The polynomial 
function generally gives a realistic fit to photographic pho- 
tometry while the exponential reproduces very well typical 
CCD photometric errors. 

The astrometry is computed using the kinematical pa- 
rameters summarized in table 4. Astrometric errors are 
assumed to follow a function of the magnitude, as the 
photometry, either polynomial or exponential. Radial ve- 
locities are computed in the same way. If needed, par- 
allaxes can also be estimated, assuming that the proper 
error function is known. 

Finally, a catalogue (or a statistics if needed) of the 
expected star distribution in any given direction of ob- 
servation is created. The output parameters for each 
simulated star include the intrinsic parameters (absolute 
magnitude, effective temperature, gravity, age, metallic- 
ity [Fe/H], U,V,W velocities computed without errors) as 
well as observational characteristics (apparent magnitude, 
colours, proper motion, radial velocity, parallax), distance 
to the sun and interstellar extinction. 

3.2. Basic applications 

The model is able to produce star counts in any given 
direction but the Galactic center (the inner bulge being 
not yet included). It produces alternatively a catalogue 
of simulated stars or statistics on any given observable or 
intrinsic property. 

The model has already been widely used (either in its 
previous version or in the present one) mainly for the fol- 
lowing purposes : 

— Predictions of star counts in different photometric 
bands from U to K in different directions to be 
compared with observational data: examples of re- 
cent comparisons betwe en model and observed star 
counts c an be fo und in Goldman et al. (2002Y'Oiha' 
J2nm^ Robin et al. (,.2000.) ,Revle fc Robin (2001a) and 
ICastellan a tOQjh in visible bands:|Krause et all ll200: " 



in the I band:^uDhv et al.' ('l996VPcrsi et al." ^1999' 



iPersi et al.„ (,2001.) and .Revle fc Robin, (,2001a) in the 
near infrared. 

Predictions in the X domain: iGuillout et al. 
have extended the predictions of the previous model 
version to ROSAT bands. It has been used to compute 
the contributions of Galactic stars to the diffuse X- 
ray ba ckground calculated for R OSAT PSPC energy 
bands ()Kuntz fc SnowdenL l200l|) . 
Predictions of kinematics and c omparis on with obser- 
vational data ( Bienavme et al.l .1992: .Chareton et al.L 



19931 lOihaet ali Il999t iRaoaport erall I2OOII: 
Soubiran et alil2003|) . 

— Predictions for use in preparing observations. The 
model allows to make simulations in different condi- 
tions, to test the limiting magnitude to be reached for 
answering a given question and to c hoose the most ef- 
ficient filters (see also the work of iLe LouarnI (|2000l) 
for studying adaptative optics possibilities with guide 
stars). 

— Determinations of e xtinction in the Galactic plane 
l|Zagurv et all Il999|) . We have also shown that the 
model can be used to map the extinction, just by play- 
ing with the extinction distribution along a line of sight 
introduced in a simulation. This will be used in the fu- 
ture to build a tri-dimensional model of the interstellar 
extinction. 

— Estimations of expected contamination of Galactic star 
in front of clusters, molecular clouds, or nearby re- 
solved g alaxies (see for e x amplelPierrg (119871) :ICananzil 
fl99i: Irhorava] et abl fl99^: IPeeg et alJ j1998^: 
'Cambrcsv ct al. (1998|): IVuongl ((200lh : lR,eikubal ((2OOII 
12002.) : Moraux (.200111 

— Simulations of proper motion selected samples to 
correc t kinematically biased selections (iRevle et al.l 

The wide variety of simulations which can be done 
allow to obtain a wide set of different type of constraints 
for the model. Hence it furnishes a wide variety of ways 
to constrain the scenario of Galactic formation. 

In order to test model hypothesis, predictions have 
been abundantly compared with star counts in many di- 
rections, photometric bands from B to K, in magnitudes, 
colours and proper motions. 

We give here just a sample of model predictions com- 
pared with observational data in magnitudes V and K 
and a colour/proper motion distribution. Figures 4 to 7 
show star count predictions in the V band as a function of 
longitude and latitude. We have added a few observational 
measurements from the literature. The agreement between 
model and data is satisfying and gives confidence for us- 
ing this model for interpolations. It should be noted here 
that the extinction has not been adjusted in the Galactic 
plane, hence predictions given at b=0° may be inadequate. 
Figures 8 to 10 show star count predictions at 2.2 /im and 
are compared with a few samples from the DENIS survey. 
In figures 5 and 8 we have not given the expected star 
counts at the Galactic center, as the inner bulge is not 
modeled. 

In figure 6 and 9, the numbers are given for longitudes 
90 and 270, for both quadrants. The counts vary slightly 
between these two directions at V>12 and |6| < 45 because 
of the warp. In the anticenter direction, as seen in figure 7, 
the effect of the flare is that the counts at latitude 10° are 
at the same level as at latitude 0. 

Figure 11 gives the contributions of various spectral 
types as a function of V magnitude in the direction 1=90° 
and b=45°, representative of what happens at medium lat- 
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itudes. Giants dominate the counts up to magnitude 11 in 
V, where main sequence stars start to be more numerous. 
The limiting magnitude between these two regimes moves 
to fainter magnitudes when the latitude goes down. 

Figure 12 shows a comparison between observed and 
modeled colour-proper motion diagram towards 1=3° and 
b=4 7° in the magnitud e range 12 < V < 18. Data are 
from lOiha et~al\ l)l996(l . The model reproduces well the 
overall distribution, the mean and dispersion in proper 
motion, and the colour distribution: the red part is dom- 
inated by the disc population, the blue one by thick disc 
and spheroid stars. 

More comparisons between model predictions and ob- 
servational data can be found in previous published pa- 
pers. Colour histograms si mulated and compa red with ob- 
servations can be found in iR obin et all ll200f]l) in 15 direc- 
tions. In iRevle et alJ (2001b) we show a velocity diagram 
made from a simulation of a p roper motion selected sam- 
ple, to be compared with the lOpoenheimer et"al 
sample of probable white dwarfs. 

3.3. Bayesian classification 

Apart from common predictions, star counts and cata- 
logues of simulated stars, the model predictions may also 
be useful for estimating Bayesian probabilities: given a set 
of observables and observational errors, by suitable simu- 
lations of the direction of interest, the model can estab- 
lish a probability classification. This classification is in- 
dependent from any spectral classification: assuming that 
the model reproduces correctly the stellar content of the 
Galaxy, it allows to perform an estimate of the probable 
class, metallicity or population of an observed star, using 
different types of observables such as magnitudes, colours, 
proper motions, radial velocities, etc. 

A pr evious use o f this Bayesian classification can be 
found in lRevle et alJ l)2002i) . In a proper motion selection 
sample of nearby dwarfs, the combined use of magnitude, 
colours and proper motions allowed to determine the prob- 
able membership of observed stars, to deduce a metallic- 
ity and a distance (as the distance estimator is metallicity 
dependent). We think that this approach can be used on 
many occasions as a complementary classification to usual 
photometric or spectral classification, in particular in the 
framework of Virtual Observatories. The limitation of this 
classification is the adequacy of the model to the data. One 
has to verify in each case that the model reproduces well 
the distribution of available observables before relying on 
such classification. 

4. Conclusion and future plans for model 
developments 

The model may suffer from misfits in fields close to the 
Galactic plane due to the crude model of extinction, which 
does not account for high or irregular absorbing clouds and 
the absence of spiral arm modelisation. While we do not 
intend to model any stellar and interstellar fluctuations 



in the Galactic plane, at least an estimate of the amount 
of extinction, and its variation along the line of sight, is 
needed to interpret observations in the disc, even in the in- 
frared. We intend to participate to the challenge of getting 
a reliable 3D model of ext inction, such as that the one of 
jPrimmel fc Spergell l)200l[) . but constrained by observed 
star counts, which, we think, are the more reliable way of 
determining the distance of the absorption features. 

At the moment, the model is well constrained by star 
counts at magnitudes between 12 and 24. At both ends, 
deviations are found and should be studied in more de- 
tail. Extrapolations at very faint magnitudes are still risky 
since the predictions strongly depend on the assumed IMF 
at low masses which are still very uncertain. Figure 4 
shows model predictions and observed star counts towards 
the Galactic pole in the V band. The overall agreement is 
very good up to V=22. At deeper magnitudes the agree- 
ment would be better by adopting a higher IMF slope for 
the spheroid (a = 2, dotted dashed fine), compared with 
the standard value of 0.5 (solid line). But we have only 
one data set to rely on at these magnitudes and we suspect 
the observations to be contaminated by galaxies at V> 22. 
The determination of the spheroid IMF at low masses re- 
lies on a few very deep star counts: space based counts are 
scarce and have poor statistics, while ground based ones 
suffer from contamination by galaxies. This point will cer- 
tainly be addressed by ongoing wide field deep surveys, 
such as the coming CFHTLS, VST or VISTA surveys, 
and from the ACS camera of the HST. 

Nevertheless the present version of the 
Besangon model of the Galaxy (available at 
http://www.obs-besancon.fr/modele/model2003.html) 
is well suited for testing Galactic structure parameters 
and evolution scenarios, and is useful for accurate pho- 
tometric and astrometric predictions in different filters 
from the visible to the near-infrared. We plan to provide 
calibrations for the SDSS system, and for the future 
GAIA photometric system, when chosen. The so urce 
counts in ROSAT bands from Guill out et"al] l|l99(i() will 
be adapted for XMM satellite bandpasses. We have also 
computed magnitudes in some ISOCAM filters which will 
be use d in the near f uture for interpreting the ISOGAL 
survey l)Omont et all l2000() . 

The mass model can be used to compute optical depths 
towards the inner Galaxy and the Magellanic Clouds, to 
be compared with constraints from microlensing experi- 
ments. This is also one of our aims to use these results to 
analyse self-consistently the measured optical depths fully 
taking into account the current knowledge on Galactic 
structure and populations. 
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Fig. 4. Star count predictions in the V band at the Galactic pole. Solid line : total with a spheroid IMF slope of 
a = 0.5, long dashed: disc, short dashed: thick disc, dotted: spheroid, dotted dashed: total with a spheroid IMF slope 
of Q = 2. Vari ous observed counts are indicated as crosses (error bars ar e the Poisson noise) from Simbad data base, 
iGilmore et all 1(1985,) . .Bok fc Basinski (,1964,) . .Chiu. (,198Q,) . .Yoshii et all l|l98l1) and Creze et al. in preparation. 
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Fig. 5. Star count predictions (stars per magnitude and per square degree) in the V band at 1=0°, for l atitudes 10° 
to 90° f rom top to bottom (20° 45° a nd 90° with solid Unes, 10°, 30°, 60° with dashed hnes). Data are from lOiha et alJ 
' l994a^ and lYoshii fc R.odgerj (jic 
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Fig. 6. Star count predictions (stars per magnitude and per square degree) in the V band at 1=90° or 270°, for latitudes 
10° to 90° from top to bottom (0°,20°,45° and 90° with sohd hnes, 10°,30°,60° with dashed hnes). For e ach latitude 
the bi furcation is due to the warp, the highest part being for 1=90, b>0 or 1=270, b<0. Data are from 'Oiha et al.| 
l|l99ti() . the DMS survey Hah et al. (1996) and in the Chandra Deep Field South from .Groeneweeen et al... (2002.) 
(longitudes and latitudes of the data are in parentheses). 
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Fig. 7. Same as figure 5 but for 1=180. Data are from lOiha erall l|l994b(l . lBo'rra k Levasl l|l986|) and lChiul IjlQSd) in 

indicated directions. 
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Fig. 8. Star count predictions in the K band (2.15 /im) at 1=0°, for latitudes 10° to 90° from top to bottom (20°, 
45° and 90° with solid lines, 10°, 30,° 60° with dashed lines). Data are from the DENIS point source catalogue in 
indicated directions. 
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Fig. 9. Star count predictions in the K band (2.15 fim) at 1=90° or 270°, for latitudes 0° to 90° from top to bottom 
(0°, 20°, 45° and 90° with solid lines, 10°, 30°, 60° with dashed hues). For each latitude the bifurcation is due to 
the warp, the highest part being for 1=90, b>0 or 1=270, b<0. Data are from the DENIS point source catalogue in 
indicated directions, and from .Kiimmel & Wagner. (,200 L*! towards the North Ecliptic Pole. 
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Fig. 11. Star count predictions (stars per magnitude and per square degree) in the V band at 1=90° and b=45° and 
contributions of the different luminosity classes. All: solid line, class II: long dashed, class III: short dashed, class IV: 
dotted, class V: dashed-dotted upper curve, white dwarfs: dashed-dotted lower curve. The giants dominate at V<11, 
main sequence stars at V> 14 in this direction. The rise of the white dwarf curve at V> 27 is due to the dark halo 
white dwarf population normalized here to 2% of the dark halo. 




